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We study ions in a nanotrap, where the electrodes are nanc 
the role of a quantum optical system which acts as a probe 
with or between nanomechanical resonators. 



Laser manipulated trapped ions are one of the prime 
examples of a quantum system, where control of coherent 
quantum dynamics, state preparation and measurement 
are achieved in the laboratory, while decoherence due to 
coupling with the environment is strongly suppressed 
These achievements are illustrated by recent progresses 
in developing ion traps for quantum computing and high 
precision measurements j2(. A key step in the future will 
be the realization and integration of ion traps with micro- 
fabricated nanostructures, such as segmented traps and 
on-chip ion traps with strong confinements] . As a new 
aspect, this opens the possibilities of developing trapped 
ions as a quantum optical system which acts as a probe 
and control, and allows entanglement with or between the 
quantum degrees of freedom of mesoscopic systems 0, Q 
- while raising interesting questions of decoherence in a 
solid state environment. 

In this paper we study ions in a mesoscopic Paul trap, 
where suspended nanomechanical resonators play the 
role of tiny trap electrodes, and act as high-Q nanome- 
chanical resonators with their own quantum degrees of 
freedom. Below we develop a model of the trapped ions 
coupled with the flexural modes of these nanomechan- 
ical electrodes. In particular, we investigate the pos- 
sibility of manipulation, preparation and measurement 
of quantum states of the flexural modes via the laser 
driven ion in the limit where the trapping frequency of 
the ion is resonant with the frequency of the nanomechan- 
ical oscillator. This setup can be generalized to ion trap 
and nanoelectrode arrays. Another application is quan- 
tum computing, where mesoscopic traps not only promise 
very strong confinement and an associated speed up of 
two-qubit quantum gates, but coupling via the nanome- 
chanical electrodes offers new ways of entangling internal 
states of ions. We also study the decoherence mechanism 
for the trapped ions due to the nanomechanical and elec- 
trical couplings which introduce quantum Brownian mo- 
tion of the electrodes and limit the ability of manipulat- 
ing the mechanical modes via the ion. 

Model for ions coupled to nanomechanical electrodes: 
We consider a schematic setup as illustrated in Fig. ^ 
where an ion is trapped between two parallel suspended 
electrodes represented by nanowires or nanotubesQ. 
When a gate voltage V(t) — Vq cos uj ac t with frequency 
LUac is applied to the electrodes via Ohmic contacts, the 
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FIG. 1: Setup. Left: ion trap with electrodes made of 
nanomechanical resonators suspended above ground and bi- 
ased with ac voltage; right: arrays of nanomechanical res- 
onators and laser manipulation of the trapped ions. 



charge on the resonator oscillates with time and leads 
to oscillating forces on the ion as well as between the 
resonators. Averaging over the fast driving frequency re- 
sults in an effectively harmonic trapping potential for the 
ion centered between the electrodes. According to Fig.Q] 
the resonators of radii ro are separated by a distance 
2do, much smaller than their length 2Lq, and located a 
distance ho (< Lq) above the ground plane. As a con- 
sequence, there will be tight trapping of the motion of 
the ion along the s-axis (see Fig. ^) which couples to 
the high-Q flexural (bending) modes of the electrodes, 
compared to looser confinement in orthogonal directions. 

In Euler-Bernoulli theoryjjj the equation of motion 
for the (small) displacement Ui(y) of the flexural modes 
of the z-th resonator (i — 1,2) is given by d 2 Ui/dt 2 + 
(EI 2 /p)d 4 Ui/dy 4 = 0, where E is Young's modulus, 
I2 the moment of inertia, and p the mass density of 
the electrodes 0. The n-th secular mode with the 
eigenfunction u in (y) has frequency uj in = EI 2 / pqf n 
with wave vector <^„. We can expand the displace- 
ment of the electrode as u.i(y) = Xi n Ui n (y) where 
Xi n when quantized as Xi n takes the role of a "po- 
sition" operator with a conjugate "momentum" opera- 
tor Pi n . The lowest flexural mode, also called the fun- 
damental mode, has the secular function itm = i(y) = 
ao [cosqn coshg^Lg — cosh qny cos qnL ) where Oq is the 
normalization coefficient and qn « (27r/Lo)0.376. Note 
the flexural modes have quadratic dependence on the 
wave vector. Other acoustic phonons at long wave length 
are the longitudinal modes and the torsional modes which 
have much higher frequencies than flexural modes. For 
nanowires with the length of the order of pin the secular 
frequencies are in the range of tens to a few hundred MHz 
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with quality factors Q > 10 4 5]. For Carbon nanotube 
0, 01 , tne Young's Modulus is E = 1 TPa; with radius 
r = 2 nm and 2L = 500 nm the fundamental mode has 
a frequency of w& = 1 GHz. 

The vibration of two suspended resonators i = 1,2 
with mass AT, and eigenmodes Wj n is described by the 
coupled parametric oscillator Hamiltonian 
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where the first term is the elastic energy and the second is 
a micromotion of the electrodes due to the electrostatic 
interaction between the oscillating charges on the res- 
onators with K mn (t) = K mn cos (to ac t). Explicit expres- 
sion for the coefficients K mn can be derived from an ex- 
pansion of the electrostatic energy for given applied volt- 
ages Vi(= V(t)), E m = -\ C lJ {{X m })V i V jl with 
Cij capacitances of the resonators as functions of the dis- 
placements Xi n . For an array of electrodes, as in Fig. ^ 
these coupling terms may (partially) compensate each 
other due to symmetry considerations. 

To derive the total Hamiltonian of the ion coupled 
to the electrodes we must consider the Coulomb in- 
teraction between the ion and the charge distribution 
on the resonators. The charge distribution on the i- 
th electrodes includes an (essentially uniform) charge 
density pi V (t) — pi V cos u ac t induced by the time de- 
pendent external voltage, and an induced charge dis- 
tribution pi q due to the presence of the ion ("im- 
age" charge). As mentioned above, we only consider 
the high frequency motion of the ion along the x- 
axis which is decoupled from the motion along y and 
z. The interaction energy between the charge of the 
ion g and the total charge density on the electrodes 
has the form J2i J ds i (Piv (*) + Piq) qo/4ne \R x e x + ye y \ 
with R x = x + (— l) l+1 do — Ui (y), where we inte- 
grate along the resonator and the denominator involves 
the distance between the ion at position xe x (compare 
Fig.GJ and the electrode with displacement m (y) e x . For 
max (| it, (y)\) <C 2Lq the assumption of a constant charge 
distribution remains a good approximation. The inter- 
action of the ion with the "image charge" is well approx- 
imated by — <7o/47re | i?^ | In ( \R X \ /r ) [SB and generates a 
small correction to the trapping and the coupling energy. 
An additional effect here is the modification of p v due 
to the small voltage of the image charge to the ground, 
which can be neglected for h ^S> d n . 

Second order expansion in x and Xij gives a harmonic 
coupling between the ion and the electrodes. The total 
Hamiltonian iJtot = H v +H± on with 7J; on the Hamiltonian 
for the ion of mass m coupled to the electrodes 
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which includes the trapped ion Hamiltonian Hamilto- 
nian with Vo(t) — 2 V ac cosuj ac t/ ln(2/i y^g + dfi/r Q do) 
modified by a (small) image charge contribution, V q — 
—ot g qo/2ireadaln(da/ro) (with a g ss 1). The last term is 
the ion-electrode coupling with gi n (t) = gi n cos ui ac t and 
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where again the image charge gives a small 
modification. Note a displacement force 

^2 gin cos (^act) X in and an electrode-electrode in- 
teraction ^2 himn cos i^act) Xi m Xi n induced by the 
static ion, which can be derived by expanding the 
Coulomb interaction, are not written in Eq. (J2J). 

We are interested in a situation where the trap fre- 
quency is near resonant with the lowest flexural mode 
n = 1, while the driving field with frequency uj ac is 
far off resonant with any of the other elastic eigen- 
modes. This, together with the fact that gi n is a 
rapidly decreasing function of n, justifies the single mode 
approximation for the electrodes: H lon — p 2 /2m + 

q Vo {t) (x - Ui(0)X p ^ /2d 2 , where X p = (aiX tl + 

a 2 A" 2 i)/2 and for symmetric electrodes and do/L a < 1, 
ol\ = ct'2 ~ 1. This Hamiltonian has the simple interpre- 
tation of a parametric oscillator where the trap center 
located at the center-of-mass (COM) displacement X p , 
providing a bilinear coupling of the ion to the electrodes. 

As a final step, we adiabatically eliminate the (micro- 
motion at the parametric drive frequency uj ac . We illus- 
trate this for the case of two symmetric electrodes with 
identical eigenmode frequencies and masses, so that the 
ion only couples to the COM mode X p . We obtain for 
the effective Hamiltonian with coupling in the rotating 
wave approximation 



Hcff = Tiio v a) a + hidbtfb + (iXb^d + h.c. 
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where b (w) are the lowering (raising) operators of the 
COM flexural mode with mass M p = 2Mi\ and a (a') 
operators for the harmonically trapped ion with fre- 
quency u„ = qoVq/ 1 \f2d\mio ac and mass m. Near res- 
onance at io v ~ u>b <C w QC , the coupling is proportional 
to the trapping frequency while decreased by the mass 
ratio m/M p : A = hu} v u\{0)^Jm/ M p . This shows the 
nanomechanical resonator with lighter mass, e.g. a sin- 
gle wall carbon nanotube, will have stronger coupling 
with the trapped ion. Typical parameters are Vq — 5 V, 
w ac /2-K — 3 GHz, d — 100 nm, h a = 10(io, r = 2nm, 
m/M p w 10" 4 , and we have w v /2-n = 1 GHz. The 
time dependent coupling K mn cos 2 uj ac t in Eq. (|TJ satis- 
fies \n mn \ < (Qc/qo) {u ac /uj )mujl < M p lo 2 at w„ ~ u b , 
where Q c — 2L$p v is the charge on the electrodes from 
the voltage source. For the above numbers, Q c w 50 
and hence the n mn term, so are the ^ g£ and h^ nn 
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FIG. 2: Energy structure of the combined system of the ion 
and the electrodes. Left: cooling circles; right:heating circles. 



terms, can be neglected due to the small mass ratio. In 
general, cantilever couplings are readily accounted for by 
transforming to a set of eigenmodes. 

The coupled oscillator Hamiltonian © allows the 
transfer of the motional states of the electrodes and 
the ion. Motional states of the ions can be manipu- 
lated by coupling internal electronic states to a laser: 
Hj = -\5a z + fie"'( a+at )(T+ + h.c, which can be added 
to Eq. ©. Here cr's are the Pauli operators of the a 
two-state system, S the laser detuning, the Rabi fre- 
quency, and 77 = ky // h/2muj l/ the Lamb-Dick parameter 
describing the laser recoil on the ion motion. Based on 
these couplings, a complete toolbox is available on the ion 
for (i) quantum state engineering, and (ii) preparation of 
pure states (ground state cooling), and (iii) state mea- 
surements using the quantum jump technique. Prepar- 
ing and analyzing motional states of the resonators are 
available via the transfer Hamiltonian when the coupling 
time 1/ A is faster than the decoherence time. 

Decoherence: Decoherence of the nanomechanical 
resonator is induced by mechanical. 7] and electrical 
noise [lfj. The mechanical noise is characterized by the 
quality factor Q due to the coupling with environments, 
such as the noise in the support of the electrodes and 
the surface of the electrodes. Electrical noisejlOf repre- 
sents the voltage fluctuations on the electrodes, including 
the noise of the resistances in the circuit, shot noise and 
low frequency noise (1/f noise). The shot noise is neg- 
ligible for a ballistic electrode, i.e. single wall Carbon 
nanotube^l- The low frequency noise is due to charge 
fluctuations in the environment and the imperfection or 
dirt on the electrodes, and is negligible at GHz frequency. 

Let Svi (t) be the noise on the electrodes defined the 
same way as Vq with i = 1,2. The voltage noise intro- 
duces: the parametric noise qo 



Sv p (t) (x — uiXp^j /2d\ 



with the spectrum Svi (u>) « $ v i ( w ) /4, and the 



linear noise qoSv m (t) [x — uiX p j /2do with <5v m — 

/3 s 2 E^H/4 and (3 g « 2\n(2h /r ) / \n(2d Q /r Q ), 
when Svi^ (t) are not correlated. Due to m/M p <C 1, 
the Xp, xX p and X p terms are a factor of m/M p or 
(rn/Mp) 2 smaller than the x 2 and x terms, so we only 
consider the x 2 and x terms. The parametric noise 
has a spectral density 71 (w) = (5x /do) 4 q 2 Sv^ (uS) /4. 



Due to the small ratio 5xo/do, 71 = 0.5Hz for the 
contact resistance at zero temperature and can be ne- 
glected. Here we have R c = h/2N±e 2 with N± the 
number of conducting channels and N±_ = 4 for Car- 
bon nanotube, and Svz (u>) — R c u)/2 at zero tempera- 
ture. The linear term has a spectral density 7 m (w) = 
(Sx Q /d Q ) g 2 5v m M/4, with Sv^ (u) = R c u/2 for un- 
corrected noise and 7 m (w„) = 0.5MHz. Note that due 
to the (6x0 /do) factor, this term is stronger than the 
parametric noise. The linear noise can be avoided when 
the two electrodes are made of one piece of metallic wire 
connecting to only one contact. Now Sv\ — Sv2, and 
7 m ~ as well as Sv^ ~ 0. Thus the condition of a de- 
coherence time longer than the coupling time of the ion 
to the resonator can be met. 

Resonator Cooling: The energy structure of the sys- 
tem is shown in Fig0 where the state \s, n, m) corre- 
sponds to the internal state at s, the motion of the ion 
at the number state n, and the vibrational mode at the 
number state m. The states \s,n,m), \s, n — 1, m + 1) 
and \s,n + X,m — l) are connected by the coupling A. 
Eliminating the internal degrees of freedom of the ions, 
the master equation in the interaction picture of the ion- 
resonator motion becomes 



dp 
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where we use the notation in [l"H 
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for the Liouvillian operators with the system operator A, 
noise spectrum 7^, and a phonon number n w . In Eq. (0} 
laser cooling of the ion contributes the £ d (a,7^,n/) 
term, where nf is the final phonon number of the ion 
in the absence any additional heating mechanisms (i.e. 
nf 0) and cooling rate 7^ = r] 2 Cl 2 /T e when r e > r)tt r 
is the rate of laser cooling. The linear voltage noise con- 
tributes the C d (a, 7 m , ns) term. The resonator damping 

and heating contributes C d (h, 7^, n^j where 7^ = LUb/Q 

and ub — ksT/hujb is the thermal phonon number of the 
supports at frequency u>f,. 

For nf = 0, the final state phonon number of the elec- 
trodes is 
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2„eff 



(r„ + 7 : ff ) (7: ff r, + 4A 2 



with r„ = uib/Q- The effective cooling rate is 7° ff = 
lu + 1m including the voltage noise with the final phonon 



number of the ion 
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n B 7m/7efi- When 7 „ 



7eff 



4 



the ion can not be cooled and hence the electrodes. As 
a necessary requirement we need j m ~ 0, i.e. the two 
electrodes in one piece (see above). In this case, when 
7° = 2A the stationary phonon number is minimized: 
(b r b) min = (k B T/Q) (4A + r„)/ (2A + T„f. With our pa- 
rameters r e = 5 MHz, n r = 300 MHz, and ^ = 2MHz, 
we have (&b) / « 0.5 with Q = 10 s and T = 4K. 

For the high trapping frequencies, the Lamb-Dicke pa- 
rameter is small with rj ~ 0.01 and as a result 7° < A. 
The ion cooling rate becomes the bottleneck for the cool- 
ing process. The problem can be avoided by trapping 
the ion at low frequency w„ « wj, ~ u> ac (and thus 
large rf). The ac driving field provides a parametric 
up-conversion of the ion to the resonator phonons, i.e. 
A — > huj v uiy/ mu}^ c l MpU>i > Lj lJ e~' lUact in Eq. ©. 

Entanglement generation: Quantum state engineer- 
ing of ion can be used to generate entangled states of 
the resonators via the coupling Eq. 0. For exam- 
ple, let us discuss the generation of | ^l? X2) + Ixi,^) 
for two resonators where and are arbitrary 
states of the electrode i = 1,2. We choose the elec- 
trodes to have different fundamental frequency ui^x 7^ 
LJb2 and coupling with the ion Aj. The initial state 
(|ti i>x) + \i,Xx)) |0i, Q2) with x for the motion of the 
ion is prepared using standard protocols of the ion trap 
qubits 12j, where the motional and the internal state of 
the ion are entangled. In a first step we tune the ion to 
the first resonator resonance lo v = LObi for a duration of 
7r/2Ai, which results in the swap \n x ,mi) — > \m x ,ni), 
and the state is now (|t, X , tpi) + \i, X , xi)) |0a). Now 
we prepare the state (It, Xx, V>i) + \l,i>x,Xi)) IO2) via a 
third internal state in the ion|l2j. Then, we tune the 
trapping frequency to uo v = for a duration of tt/2X 2 , 
which results in |n 2 ,,m 2 ) — > \m x ,n2), and the state is 
(\1,0 x ,ipi,X2) + \l,0x,Xl,ip2))- Now we rotate the inter- 
nal state by a ir/2 pulse: |f) -> |t + I) and |J.) -> |f - 1); 
then detect the internal state. The detection generates 
states \ipi, X2) ± 1x1,^2) depending on the detected in- 
ternal states; and entanglement is transferred from the 
ion to the resonators. This method can be applied to 
generate for example |0i, n 2 ) + |ni, O2) as entangled state 
of the two electrode modes. The entanglement can be 
achieved on a time scale of 7r/Ai,2 ~ 50nsec. Given a 
decoherence time of 1MHz, a fidelity exceed 0.9 can be 
achieved. 

Quantum Computing: Mesoscopic traps promise trap 
frequencies significantly higher than those of present ion 
traps and an associated speed up of two-qubit gates in 
ion trap quantum computing (Fig. The standard 2- 
qubit protocols based on, for example, entanglement via 
a phonon data bus or a push gate are readily adapted 
to the present case. In the second case, gate times of 
the order of nanoseconds seem possible for the numbers 
above. A second possibility is entanglement via exchange 
of phonons of one of the collective ion-electrode modes 
of the system. Additional noise is, however introduced 



by the mechanical motion. The decoherence of the ion 
motion when in resonance with one mode lo v = loh is 
T res — ksT/Q which is 1MHz with Q = 10 5 and assum- 
ing the same quality factor for two electrodes. When off 
resonance \u>n — lo w \ > ksT/Q, the decoherence rate is 



E 
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so that the decoherence rate is T m m 
(8ulm/9M p ) (k B T/Q) - 2n- 100 Hz with u ix = 500 MHz 
and u„ = 1 GHz. The small mass ratio m/M p and the 
off resonance protect the ion from decoherence. The 
decoherence due to the electrical noise is, dominated 
by 471 (0) = 2-7T • 0.1Hz assuming a low frequency noise 
(qoSvf (0)) = In sec. Hence, the decoherence time 
is much longer in the nano trap. This shows that the 
smallest charge box - single ion - coupling with the 
nanomechanical resonators not only presents an effective 
knob for quantum features of the flexural modes, but 
also increases the speed of ion trap quantum computing 
by several orders of the magnitudes. 

Acknowledgments: We thank M.S. Dresselhaus, B.I. 
Halperin, D. Leibfried, L.S. Levitov, M.D. Lukin, A.S. 
S0rensen, and W. Zwerger for helpful discussions. This 
work is supported by the Austrian Science Foundation, 
European Networks and the Institute for Quantum In- 
formation. 



[9 
[10 

[11 
[12 



J. I. Cirac and P. Zoller, Phys. Today March vol.57, 38 
(2004). 

M. D. Barrett et al, Nature 429, 737 (2004); M. Riebe 
et al, Nature 429, 734 (2004) 

D. Kielpinski et al., Nature 417, 709 (2002); J. I. Cirac 
and P. Zoller, Nature 404, 579 (2000). 
D.J. Wineland et al, J. Res. Natl. Inst. Stand. Technol. 
103, 259 (1998), see sect. 6.4, pp319; L. Tian and et al, 
Phys. Rev. Lett. 92, (2004) 247902; A. S. S0rensen et al, 
Phys. Rev. Lett. 92, 063601 (2004). 

M.D. LaHaye et al, Science 304, 74 (2004); R.G. Knobel 
and A.N. Cleland, Nature 424, 291 (2003). 
A. Husain et al, Appl. Phys. Lett. 83, 1240 (2003); 
P. Poncharal et al, Science 283, 1513 (1999); B. Babic 
et al, Nano Lett. 3, 1577 (2003). 

A. N. Cleland and M.L. Roukes, J. Appl. Phys. 92, 2758 
(2002). 

R. Saito et al. , Physical Properties of Carbon Nanotubes, 
Imperial College Press, London (1998); MRS Bulletin 29, 
Advances in Carbon Nanotubes, April 2004. 

B. E. Granger et al, Phys. Rev. Lett. 89, 135506 (2002). 
Ya. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 
(2000). 

J. I. Cirac et al, Phys. Rev. A 46, 2668 (1992). 

C. K. Law and J. H. Eberly, Phy. Rev. Lett. 76, 1055 
(1996); J. I. Cirac and P. Zoller, Phys. Rev. Lett. 74, 
4091 (1995). 



